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(D Abstract 

Q We visualize the fundamental property of pQCD: the smaller size of the colorless quark-gluon 

O configurations leads to a more rapid increase of its interaction with energy. Within the frame of 

^ the dipole model we use the kt factorization theorem to generalize the DGLAP approximation 

and/or leading ln(,To/x) approximation and evaluate the interaction of the quark dipole with a 
target. In the limit of fixed and a; ^ we find the increase with energy of transverse momenta 
of quark(antiquark) within the qq pair produced by the strongly virtual photon. The average 
is evaluated aualyiically within the double logarithmic approximation. We demonstrate that the 
r^. invariant mass^ of the qq pair increases with the energy as Mq{xo/x)^, where A ~ OAagNc/iT for 

^^^^ transverse photons, and as ~ Mq exp0.17[(4a5A^c/7r) log(xo/x)]^/^ for longitudinal photons, where 

> PS 0. 7(^2 at the energies of the order sq ~ 10^ GeV^ (xq ~ 10 ^. The magnitude of the effect 

depends strongly on the small x behavior of the gluon distribution. Similar pattern of the energy 
^ dependence of is found in the LO DGLAP approximation generalized to account for the kt 

CO factorization. We discuss the impact of the found phenomenon on the dependence of the coherence 

^-H length on the initial energy and demonstrate that the shape of the final hadron state in DIS has 

the biconcave form instead of the pancake. Some implications of the found phenomena for the 
O hard processes in pp collisions are discussed. 
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I. INTRODUCTION. 



A dipole model developed in ref. [T] , cf.also |21 El SI El E] is the generalization of the parton 
model to the target rest frame description. It accounts for the effects of the and ln(a;o/a;) 
evolutions. It also provides the solution of the equations of QCD in the kinematics of fixed 
and not too small x = Q^/iy but oo. The characteristic feature of this solution is 

the approximate Bjorken scaling for the structure functions of DIS, i.e. the two dimensional 
conformal invariance for the moments of the structure functions. In this approximation as 
well as within the leading ln(xo/a;) approximation, the transverse momenta of quarks within 
the dipole produced by the local electroweak current are restricted by the virtuality of the 
external field: 

>^lcD<Pt<iQ')/^- (1-1) 

The aim of the present paper is to demonstrate that the transverse momenta of (anti) quark 
of the qq pair produced by a local current also increase with energy and become larger than 

at sufficiently large energies. Technically this effect follows from the more rapid increase 
with the energy of the pQCD interaction for smaller dipole and kt factorization. Let us note 
that this phenomena is different from the well known Lipatov diffusion. The latter means 
that within the leading log asln{xo/x) approximation the parton transverse momenta are 
increasing with energy in the center of rapidity as ln^(p^/p^Q) oc ln(s/so) as a result of the 
diffusion in the space of transverse momenta[7]. 

Within the double logarithmic approximation we evaluate analytically both the maximum 
in the distribution over the invariant masses of the qq pair which contribute to the transverse 
and logitudinal total cross section of DIS, and the corresponding average transverse momenta 
squared. 

Consider first the case of the longitudinal photons. Then the position of the maximum 
increases with energy as 

= MIl{s/sqY^^^^'^'^'\ (1.2) 
Eq, 1.2 is derived in the approximation -C <^ s which is self-consistent at sufficiently 



high energies. One can see from this expression that that transverse momenta of quarks 
increase with the energy since = p1/z{l — z), and the configurations with 2; = 1/2 
dominate at sufficiently high energies. 

The dependence of the average quark transverse momenta on energy is calculated below 
numerically within the double logarithmic approximation and/or within the LO DGLAP 
and BFKL approximations. For certainty we define average transverse momentum of quark 
as corresponding to the median of integral for the total cross section. Within the double 
logarithmic approximation to the cross section initiated by longitudinal photon we obtain: 

Ml ~ Apl ~ 0.7Q'exp(0.17((a,Ar,/7r)log(a;o/x))°-^^). (1.3) 

Here xq ~ 0.01. The analysis was done in the interval s = 10"^ to s = 10^^ Gev^. Note that 
the derived rate of the increase with the energy of the characteristic scale does not depend 
on the external virtuality Q^. However, Mq depends on the normalisation point in xq and 
Qq. It is worth emphasizing that since we are interested here in the proof of the rise of 
the transverse momenta in the current fragmentation region, we carry for the illustration, 
the calculations over a very wide spectrum of energies s ~ 10"^ 10^^ GeV^. The detailed 
calculations for the realistic energies have been carried in the LO approximation using the 
CTEQ5L gluon pdfs [9l [in]. Quahtatively they produce similar results although depended 
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on chosen extrapolation to small x. In particular the CTEQ6L parametrization leads to a 
significant suppression of the effects discussed in the paper. 

Similar results were obtained for the transverse photons. In this case we were able to 
carry out an analytical calculation for the invariant mass distribution maximum for the 
symmetric configurations and found that it rapidly increases with energy: 

~ Ml{slso)T^^''^'^^'\ (1.4) 

The analytical results has been obtained in the kinematics: <^ ^ s. It is well 
known however that in the case of the transverse photons a major role in a wide kinematical 
region is played by qq configurations where one of the partons carries most of the plus 
component of the photon momentum. With increase of the energy the role of asymmetric 
configurations is reduced since their contribution grows with energy more slowly. In order 
to take into account the asymmetric configurations we have made a numerical calculation 
of a transverse cross-section in the interval s = 10^ — 10^^ GeV^, and obtained: 

M| ~ 0.7g2(a;o/x)°-^"^^^/", (1.5) 

xo ~ 0.01. 

If we take into account the increase of the transverse momenta of the dipole p'f with 
energy within the framework of the dipole model and the kt factorization theorem we are 
lead to the generalization of the DGLAP [8J and BFKL [7] approximations which is done in 
the paper within the LO approximation. 

The rapid increase of the characteristic transverse scales in the fragmentation region 
has been found first in ref. |TT1 [121 [131 E], within the black disk (BD) regime. Our 
new result is the prediction of the increase with energy of the jet transverse momenta 
in the fragmentation region, in the kinematical domain where methods of pQCD are still 
applicable. This effect could be considered as a precursor of BD regime indicating the 
possibility of smooth matching between two regimes. 

As the application of obtained results we obtain that in pQCD 

adx, Q')/aT{x, Q') oc (gV4p?) oc {Q^/s)\ (1.6) 

i.e. this ratio should decrease as the power of energy instead of being 0{as). 

The increase of the parton momenta in the DIS in the current fragmentation region leads 
to the change of many characteristics of high energy processes. We find that the coherence 
length of the DIS processes increases with energy within pQCD as 

oc {l/2mM){s/Q^f-\ (1.7) 

i.e. slower than in the parton model (l/2mjva; - the loffe length). This is the because the 
coherence length for a given process follows from uncertainty principle: 

h = {sl2mN)l{M\s) + Q^), (1.8) 

where M^(s) oc p1{s) is the typical important in the wave function of photon in the 
target rest frame and pt is the transverse momentum of constituents in the wave function of 
photon. This result has the implication for the space structure of the wave packet describing 
a rapid hadron. In the classical multiperipheral picture of Gribov a hadron has a shape of 
a pancake of the longitudinal size l//x (where is the soft scale) which does not depend on 



3 



the incident energy fT5] . On the contrary, we find the biconcave shape for the rapid hadron 
in pQCD with the minimal longitudinal length for small impact parameter b decreasing with 
increase of energy and being smaller for nuclei than for the nucleons. 

The paper is organized in the following way. In section 2 using the technique first in- 
troduced in QED by V.Gribov [12], we rewrite the formulae of the dipole model for the 
inelastic cross-section of DIS in the form of the spectral representation over invariant masses 
for both longitudinal and transverse photons, kt factorization [TTl [18] is explicitly fulfilled 
in this representation. 

The analysis of these formulae predicts increase with energy of transverse quark momenta 
in the current fragmentation region. In section 3 we use the double logarithmic approxi- 
mation for the amplitude for the interaction of quark dipole with the target, to evaluate 
the increase with the energy of the quark transverse momenta in the current fragmentation 
region. In section 4 we study the dependence of coherence length on the collision energy. In 
the section 5 we explain that in pQCD rapid hadrons and nuclei look like bi-concave lenses. 
Finally, in section 6 we discuss the possible applications of our results to pp, pA collisions 
at the LHC. 



II. THE TARGET REST FRAME DESCRIPTION. 



Within the LO approximation the QCD factorization theorem allows to calculate the 
total cross section of the longitudinally polarized strongly virtual photon scattering off a 
hadron target through the convolution of the virtual photon wave function calculated in 
the dipole approximation and the cross section of the dipole scattering off a hadron. In the 
target rest frame the cross section for the scattering of longitudinally polarized photon has 
the form [H El [HI [20] : 



cr(72 + T^X) 



127r2 



(2.1) 



(2.2) 



Here a is the dipole cross-section operator: 

a = (47rV3)a.(Pt')(-A) • xG{x = (M^ + Q^)/s, M^)), 

A is the two dimensional Laplace operator in the space of the transverse momenta, and 

M' = ip',+ml)/zil-z), (2.3) 

is the invaiant mass squared of the dipole. In the coordinate representation a is just a mul- 
tiplication, but not a differential operator. In the leading ln(a;o/a;) approximation a similar 
equation arises where the cross section is expressed in terms of convolution of impact factor 
and unintegrated gluon density. In practice, both equations should give close results. Using 



an integration over parts over pt it is easy to rewrite Eq. |2.1| within the LO accuracy in the 
form where integrand will be explicitly positive: 



a(7l + T^X) 



where 



/a.A..(V.,(..)|/(..rf,|v.,(.,, 

/ = (47rV3)a,(p?)a;G(5,M2). 



(2.4) 
(2.5) 
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In the derivation we use boundary conditions that photon wave function is neghgible at 
pf ^ oo and that the contribution of small pt is the higher twist effect. 

The cross section of the interaction of the longitudinal photon can be rewritten in the 
form of spectral representation by explicitly differentiating the photon wave function : 



= Svr^ — ^--^ ^ / dM'a,{My4)---f---- ■ g{x,M'). (2.6) 



— — / dM^as(M^ 4:)- — 

12 J ' ^(M2 + g2)4 

Here = 4/3 for the colorless dipoles build of color triplet constituents, and = 9/4 for 
the gluonic dipoles. 

The spectral representation of the electro-production amplitude over is a general 
property of a quantum field theory at large energies where the coherence length significantly 
exceeds the radius of the target T |2Tl[22]. The pQCD guarantees additional general prop- 
erty: the smaller size of the configuration in the wave function of projectile photon leads to 
the smaller interaction with the target but this interaction more rapidly increases with the 
energy. In the NLO approximation the structure of formulae should be the same except the 
appearance of the additional qqg, ... components in the wave function of photon due to the 
necessity to take into account the QCD evolution of the photon wave function [,19j. 

The similar derivation can be made for the scattering of the spatially small transverse 
photon. In this case the contribution of small pt region (Aligned Jet Model contribution) is 
comparable to the pQCD one. To suppress AJM contribution we restrict ourselves in the 
paper by the region of large pf and sufficiently small x where pQCD contribution dominates 
because of the rapid increase of the gluon distribution with the decrease of x. 

The pQCD contribution into the total cross-section initiated by the transverse photon 
has the form: 
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+ (1-^)2 (M^ + g^ 



(2.7) 



Here while doing the actual calculations we introduced a cut-off in the space of transverse 
of transverse momenta M'^z{l — z) >u,u^ 0.3 GeV^. 



III. THE DOUBLE LOGARITMIC APPROXIMATION. 

In this section we analyze the new properties of the pQCD regime within the double log 
approximation. The advantage of this approximation is that it will enable us to perform 
some of the calculations analytically/ semianalytically. We will present the detailed numerical 
results for the pattern observed in the complete LO/NLO in a later more detailed publication 
where we will demonstrate that qualitative though not numerically the results are the same 
(see however some preliminary results in LO below). 

We find that at sufficiently large energies the characteristic invariant mass of the system of 
constituents produced by the electromagnetic current is not but much larger-M2(s, Q'^). 
Thus in the calculations of of the high energy processes the effective virtuality is but 
not Q^. 



5 



In the double logarithmic approximation the structure functions are given by [23] 

g{x,Q^) = f dj/{2m){x/xoy-\Q'/Qly^^\ (3.1) 



where the anomalous dimension is 

<J - 1) 

To simplify the calculation we assume, the initial condition for the evolution with Q^: 

g{x,Ql) = 5{x-l). (3.2) 
In the saddle point approximation one finds |23j : 

9{x, Q') = i^g(;^/;)3/4 \/4«s(gg)iV,/vr log(gVgg) log(xo/x) (3.3) 

Structure function of a hadron is given by the convolution of this kernel with nonperturbative 
structure function of hadron in the normalization point = Qo- 

In the following we shall neglect the pre-exponential factor, since absolute value of g as 
well as the pre-exponential factor weakly influence the transverse scale, and its evolution 
with energy: 

gix, Q2) = exp ^AasiQl)Nj7T\ogiQyQl)\ogixo/x). (3.4) 

A. Energy dependence of the quark transverse momenta for fragmentation pro- 
cesses initiated by longitudinal photon. 

We shall find analytically the scale of the transverse momenta in the limit when s 3> 
^ . For certainty we restrict ourselves to the contribution of light quarks. 

At large the cross-section for the scattering of the longitudinal photon is dominated 
by the contribution of the spatially small dipoles, so it is legitimate to neglect the quark 
masses. In this limit the cross section is proportional to 

aiocQ'^ J dM^n{M\s,Q^), (3.5) 



where the function n{M'^, s, Q^) is given by Eqs. 2.6, 2.7 
n(M2,.,Q2) = a,(MV4) 



(M2 + Q2) 



X exp(v^4a,(g2)(iV,/7r)ln(MVM2)ln((V(M2 + Q2)((M2 + g2)/so)), 

(3.6) 

Here we keep only large terms depending on (we do not write here explicitly the 
independent overall normalization factor irrelevant for the calculations below). 
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Let us show that the maximum of n{M'^, s, Q^) increases with the energy. At very high 
energies n is proportional to 

n ~ exp(lna,(MV4) + hgiM^Q^) - 41og((Q2 + M^)/Q^) 



+ v'4a.(iVc/vr)(ln(s/so) - log((Q2 + M^)/{Q^ + M^)) ln(MVM2)). 

(3.7) 

In the hmit of fixed but very large energies, log(s/so) > log(((5^ + + Mg )). Let 

us assume that for the maximum: ^ Q^. We can find the maximum of the expression 



3.7 under this assumption analytically, and then check that this assumption is indeed self- 
consistent. As a result we can rewrite Eq. |3.7| as 



n ~ exp(ln(a,(MV4) - 3 log(MVQ') + ^j^a,{Ql){NJ^x){\Yl{s/so)HMyMi)). (3.8) 

Differentiating the argument of the exponent over log(M^/MQ) we obtain the equation for 
the maximum: 



l/ln(MV4M2) + 3 = {l/2){l/ \n{M^/Ml)^Aas{Ql){NjT^) \n{s / s^) / HM^ M^) . (3.9) 

Neglecting the small first term we find: 

= M2(s/so)"=(^^/"^/^ (3.10) 

Here Mq ~ and sq ~ Q^. We will refer to this extremum value of as . 
At the extremum n oc (a/(Mf /4)/Mf)exp((A^c/7r)(as/3) ln(s/so)). Therefore 

^M=M? ^ as{Ml/A){Q'/Mt){eM{Nc/r^){as{Ql)/^)Hs/so)) (3.11) 

However, the position of the maximum of the integrand is not sufficient to characterize the 
relevant transverse scales as a large range of is important in the integrand. (Calculation 
of second derivative shows that dispersion over is large. The width of the distribution over 
ln(M7M2) is ■y(2/3)log(MVM2).) Hence we need to determine M range which gives most 
of the integrand support. For certainty, we define the range of < which provides 
a fixed, say, 50% fraction of the total perturbative cross-section. Let us estimate how this 
scale increases with the energy in the double log approximation. First, let us consider the 
total cross section. The upper limit u of integration over is determined by the condition 
< s. 

For certainty we choose upper limit of integration as 

< M\^, = 0.2s, (3.12) 

although the result of numerical calculations is insensitive to the upper bound because 
essential are significantly smaller. 

Let us first calculate the median scale semianalytically. Within the double logarithmic 
approximation, and assuming that the conditions log(s/so) ^ log((Q^ + M^)/(Q^ + Mq), 
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is still valid for the relevant M^, the integral for the cross section can be written similar to 



Eq,3.8 



r\og{u/M^) 

a{u) = {Q^/M^) / ciln(MVMo2)a,(MV4)exp(-21n(MVM2) 







(3.13) 



Here u is the upper cut-off in the invariant masses. Introducing the new variable t = 
log(MVMo2), we obtain: 

a{u) = (Q^M^) ^"^"^ dtas{tM^/4) exp(-2t + ^{4asNjn)log{s/so)t), (3.14) 



where k{u) 
lar to Eq. 



3.14 



n{u/sQ). The integral for the total cross-section is given by the equation simi- 



with the upper integration limit being replaced by k{s) = y ln(0.2s/so). The 
integral 3.14 is actually the error function j23|, which can be easily evaluated numerically. 
Requiring that it gives one half of the cross section we find 

~ Mo2(s/so)°-'^"=^=/". (3.15) 

Evidently, for sufficiently large s our initial assumption log(M^/MQ) S> log((5^/Qo) 
fully self-consistent. This is because the decrease of n with due to terms in the 



integrand of Eq. |2.6| is partially compensated by the rising exponential, giving a relatively 
slow decrease of n to the right of its maximum. 

Note that the rate of the increase of with s is much higher than for Mf due to the 
slow decrease of the integrand with M^. The cross section of jet production with at this 
interval also increases with the energy as 

^^Im2=m2 ~ [s/so) ^' . (3.16) 

In order to understand the dependence of the median scale on both the energy and 
quantatively we made numerical calculation of the characteristic transverse momenta using 
the DGLAP double log structure function. We find that the increase rate of the transverse 
momenta indeed does not depend on the external virtuality Q^. Considering the wide interval 
of energies and s = 10^ -r- 10^^ GeV^, and 20 < < 200 GeV^ we obtain the approximate 
formulae: 

~ 0.7g2) exp(0.17((4a,iVe/vr) log(xo/x))°-^^). (3.17) 

We give this estimate only for illustrative purposes, since the double logarithmic approxi- 
mation is semirealistic only. Still our results indicate that for external virtualities < 100 
GeV^ and energies which can be reached at LHeC the onset of a new pQCD regime may 
take place. 

The cross section of the jet production at this scale also increases with the energy as 
da . 1 /(x) 



^M^\M^=M? gn + 0.7/(a;) 
where 

fix) = exp(0.17((4a,iV,/7r) log(a;o/x))°-^^). (3.19) 



-G{x{l + 0.7f{x),0.7Q'f{x), (3.18) 
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B. Transverse photon: the chjiracteristic transverse scale in the photon fragmen- 
tation region. 



It is well known that the main difference between the longitudinal and transverse structure 
functions in the DIS is the presence of the strongly asymmetrical in z configurations due 
to the presence of the {^z{\ — z))'^ multiplier in the spectral density. As a result there is a 
competition between two effects. One is a slower decrease of the spectral function with Af^ 
(by the factor M^/Q^), leading to the more rapid increase of the characteristic transverse 
momenta for the symmetric configurations. The second effect is presence of the asymmetric 
[z — > 0) configurations which are characterized by the small transverse momenta hi for a 
given invariant mass M^. For such configurations the rate of increase of the gluon structure 
function with energy is small. 

Let us first show that the transverse momenta increase rapidly for symmetric configura- 
tions. The spectral representation for the transverse photon for symmetric configurations is 
proportional to 



(M2 + g2)4 

X ^4a,(7V,/7r)(ln(s/so) - log((Q2 + M2)/(Q2 + Ml)) \n{M^/Mi)). 

(3.20) 

Similar to the case of the longitudinal photon we obtain for high energies, when ^ Q^, 
the dependence of the maximum of n on energy: 

Ml ~ Mo'(s/so))"'^'^'=/''^/^ (3.21) 

i.e. the increase rate is twice as fast as compared to the case of longitudinal photon. Mf 
increases with s at high energies and thus the condition ^ is perfectly self-consistent 
at very high energies. 

The jet cross section at the maximum of the curve also increases as 

^|m^=m?~(^Ao)"^^=/^'^ (3.22) 

In addition we calculate the total cross-section in the same approximation semi- 
analytically getting the error function and obtain the rate of increase (s/so))'^'^'^^^"*'^°^'^\ 
which is twice that for the longitudinal case. 

However, as we mentioned above, a considerable contribution of the nonsymmetrical 
configurations has the opposite effect. In order to take these configurations into account we 
performed a numerical calculation using the gluon distribution function within the double 
log approximation. The result is that the characteristic median scale increase like 

M2(s/so)°-'(^"^^=/'^) (3.23) 

The value of the exponent is 0.12 for the beginning of the studied energy range s ~ lO^'-hlO^^ 
GeV^, and decreases to 0.09 at the upper end (for typical = 0.25. Thus the rate of the 
increase with the energy is approximately the same as for longitudinal photons for not very 
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high energies. For very high energies the symmetric configurations win over asymmetric 
ones, leading to twice as rapid increase of the transverse momenta than in the longitudinal 
case. 

The precise determination of the scale Mq(Q^) is beyond the accuracy of this paper. 
Effectively we obtain the dependence ~ 0.7Q2(xo/x)°-^(^°''^'=/''). 

One can also estimate the rate of the increase of the jet production cross section : 

h{x) = (a;o/a;)°-^(^"»^^/") (3.25) 

We found a rapid increase of the jets multiplicity. Thus the rate of the increase with energy 
of the transverse momenta of quarks in the current fragmentation region for transversely 
polarized photon is significantly more rapid. Consequently we find that ctl/o't ~ Q'^/M'^ 
being numerically small should slowly decrease with energy at sufficiently high energies . 

We conclude that it is possible to show analytically that for very high (asymptotic) 
energies the relevant invariant masses extend well beyond and increase with the energy. 

The direct numerical calculation of the scale shows that the rate of increase is inde- 
pendent of external virtuality. 

C. The leading logarithmic approximation. 

The above results were obtained in the double logarithmic approximation. It is also 
possible to carry out the numerical calculation in LO approximation using the CTEQ5L 
gluon distribution functions |25]. In this approximation the median scale still increases 
as ~ 0.7Q'^{xo/x)^, where xq ~ 10~^, and A ~ 0.06 for longitudinal and A ~ 0.08 
for transverse photons. Although this increase is quite slow, the rise of momenta is not 
negligible: for energy increase from 10^ to 10^ GeV^ the scale increases by a factor ~ 1.5. 
The use of CTEQ6L will decrease the considered effects. 

IV. THE COHERENCE LENGTH. 

In the previous sections we determined the energy dependency of the effective transverse 
scale at high energies which allowed us to evaluate coherence length. The coherence length 
Ic corresponds to the life-time of the dipole fluctuation at a given energy in the rest frame 
of the target. The original suggestion of the existence of the coherence length in the deep 
inelastic scattering was flrst made by loffe, Gribov and Pomeranchuk |2T1 [26]. It was found 
already in the sixties within the parton model approximation by loffe [27] that the coherence 
length at moderate xb is Ic ~ l/2mjvXB i.e. it linearly increases with energies. In pQCD 
coherence length 

L,= (l/2m^x)(so/s)^ (4.1) 

Less rapid increase of Lc with energy has been found before in the numerical calculations of 
structure functions in the target rest frame accounting for evolution of structure functions 
. [291EQ] 

It is worth noting that the discussed pattern of the energy dependence of the coherence 
length leads to a change of the structure of the fast hadron wave function as compared 
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to the Gribov picture [15] where the longitudinal size of the hadron is determined by the 
wee parton cloud and energy independent ~ l//x. Here /i ~ .3 ^AGeVjc is the soft 
mass scale. On the other hand a slower than Xjm^x rate of the increase of the coherent 
length with energy leads to a decrease of the longitudinal size of the hadron with energy. 
The typical size is determined by the BD momentum at a given impact parameter for a 
particular energy. Moreover since the BD momentum is larger for small impact parameters 
the nucleon has a form of a double concave lens. It is of interest also that for the zero impact 
parameter the longitudinal size of a heavy nucleus is smaller than for a nucleon. 

V. THE FORM OF NUCLEON, NUCLEUS IN DIS 

Our results have the important consequences for the transverse structure of the hadrons 
and nuclei. 

Let us consider the longitudinal distribution of the partons in a fast hadron. As it was 
already mentioned in the previous section, in the parton model the longitudinal spread of 
gluonic cloud is ~ l//i for the wee partons and is much smaller than for harder partons, 
with Lj, ~ 1/xPh for partons carrying a finite x fraction of the hadron momentum [TJ]. The 
picture is changed qualitatively in the limit of very high energies when interactions reach 
BD regime for kt 3> In this case the smallest possible characteristic momenta in the 
frame where hadron is fast are of the order kt{BDR) which is a function of both initial 
energy and transverse coordinate, h of the hadron. Correspondingly, the longitudinal size is 
~ l/kt{BDR) -C l//i. Since the gluon parton density decreases with the increase of h the 
longitudinal size of the hadron is larger for large 6, so a hadron has a shape of biconcave 
lens, see Figs. 1, 2. 

We depict the typical transverse structure of the fast nucleon in Fig. 1. We see that it 
is drastically different from the naive picture of a fast moving nucleon as a flat narrow disk 
with small constant thickness. 

Consider now the case of the DIS on the nuclei. 

Consider first the case of external virtualities of the order of several GeV. In this case the 
shadowing effects mostly cancel the A^/'^ for a given impact parameter, h [19] and the gluon 
density in the nuclei is comparable to that in a single nucleon for 5 ~ 0. Consequently over 
the large range of the impact parameters the nucleus longitudinal size is approximately the 
same as in the nucleon at 6 ~ 0. 

Consider now the case of the large external virtualities > 40 GeV^. In this case the 
leading twist shadowing is small, and the corresponding gluon density unintegrated over b 
is given by a product of a nucleon gluon density and the nuclear profile function: 

T(6)= jdzp{h,z), (5.1) 

where the nuclear three-dimensional density is normalized to A. We use standard Fermi step 
parametrization [28] 

pir) = C(A) — — , R = l.lA^/^fermi, a = 0.56fermi. (5.2) 

^ ^ ^ ^1 + exp((r - R)/a ^ ' 

Here r = ^/z'^ + b"^, and A is the atomic number. C(A) is a normalization factor, that can 
be calculated numerically from the condition / d^rp{r) = A. At the zero impact parameter 
T{b) ^ 0.5^1/3 for large A. 
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The dependence of the thickness of a fast nucleus as a function of the transverse size is 
depicted in Fig. 2 for a typical high energy s = 10^ GeV^, = 40 GeV^. We see that the 
nuclei also has a form of a biconcave lens instead of a fiat disk. 

Note that this picture is very counterintuitive. We see that the nuclei is thinner than a 
single nucleon, i.e. the thikness of a nuclei is smaller than for a nucleon although we have 
~ A^^^ nucleons at the central impact parameter. In other words the longitudinal extend of 
a system is smaller than the the longitudinal extend of its constituents. The explanation of 
this BDR phenomenon is straightforward. For a given impact parameter b, the longitudinal 
size of a heavy nucleus l/k\^\BDR) < l/k[^\BDR) since the gluon distribution function 
in the nuclei (^^(x,^) > GN{x,b). So a naive classical picture of a system build of the 
constituents being larger than each of the constituents is grossly violated. This situation 
is in some respects analogous to the phenomenon of color transparency/existence of the 
point-like configurations). Thus we have a paradox: in the fast reference frame the nuclei is 
much thinner than any of its constituents. 




FIG. 1: 3D image of the fast nucleon at s = 10^ GeV^ and = 40Ge\'^. 

The resolution of the paradox in the BD regime is quite simple: the soft fields of individual 
nucleons destructively interfere cancelling each other. So a naive classical picture of a system 
build of the constituents being larger than each of the constituents is grossly violated. 

VI. EXPERIMENTAL CONSEQUENCES. 

The current calculations of cross-sections of hard processes at the LHC are based on the 
use of the DGLAP parton distributions and the application of the factorization theorem. 
Our results imply that the further analysis is needed to define the kinematic regions where 
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FIG. 2: 3D image of the fast heavy nucleus (gold) at s = 10^ GeV^ and = 40GeV^ 

one can use DGLAP distributions. We showed in the paper that for DIS at high energies 
there are kinematic regions where one is forced to use a kt factorization and the dipole model 
instead of the direct use of DGLAP. A similar analysis must be made for the pp collisions at 
LHC The expected effect is the increase with energy of the probability of the small dipoles 
in the wave function of proton [32] . Quantitative analysis of this problem will be presented 
elsewhere. 

The hard processes initiated by the real photon can be directly observed in the ultraphe- 
ripheral collisions [31]. The processes where a real photon scatters on a target, and creates 
two jets with an invariant mass M^, can be analyzed in the dipole model by formally putting 

= 0, while is an invariant mass of the jets. In this case with a good accuracy the 
spectral density discussed above will give the spectrum of jets in the fragmentation region. 
Our results show that the jet distribution over the transverse momenta will be broad with 
the maximum moving towards larger transverse momenta with increase of the energy and 
centrality of the 7^4 collision. 

Finally, our results can be checked directly, if and when the LHeC facility will be built at 
CERN. One of us B.Blok thanks S.Brodsky for the useful discussions of the results obtained 
in the paper. This work was supported in part by the US DOE Contract Number DE- 
FG02-93ER40771 and BSF. 
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